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ROLE  OF  THE  SECTIONALIZED  FOURIER  TRANSFORM 
IN  HIGH-SPEED  COHERENCE  PROCESSING 


INTRODUCTION 

With  the  advent  of  the  Fast  Fourier  Transform  (FFT)  algorithm  in  the  mid-1960's  [1-3]  and  the 
corresponding  advances  in  digital  computer  architecture  (in  particular,  in  array  processors),  giant  strides 
have  been  made  in  the  rapid  computation  of  complex  functions  which  earlier  were  considered  impracti¬ 
cal.  More  recently,  interest  has  evolved  in  techniques  for  estimating  the  magnitude-squared  coherence 
(MSC)  function  [4-16],  and  in  means  for  rapidly  mapping  this  estimate  over  the  two-dimensional  ambi¬ 
guity  plane  [17-20].  These  techniques  invariably  involve  computing  the  Fourier  Transform  of  the 
relevant  temporal  functions  in  a  piecewise  or  sectionalized  manner,  and  algorithms  for  computing  MSC 
estimates  (using  the  FFT  and  modem  array  processors)  have  proven  highly  successful.  Unfortunately, 
knowledge  of  these  techniques  (in  their  entirety)  has  not  been  widely  disseminated.  Nor  has  the  role 
of  the  sectionalized  Fourier  Transform  in  signal  processing  applications  been  thoroughly  understood  by 
the  user  community.  This  report  is  therefore  devoted  to  developing  the  fundamental  role  that  the  sec¬ 
tionalized  Fourier  Transform  plays  in  temporal  signal  processing,  and  to  developing  a  high-speed  algo¬ 
rithm  for  estimating  the  normalized  correlation  envelope  (NCE)  function  over  the  two-dimensional 
ambiguity  surface. 


SIGNAL  TRANSFORMATION 

Two  approaches  will  be  taken  to  demonstrate  the  role  of  the  sectionalized  Fourier  Transform 
(SFT)  in  temporal  signal  processing.  In  either  approach  it  is  shown  that  the  SFT  can  serve  to  simply 
baseband,  filter,  and  sample  a  narrowband  signal.  Due  to  the  filtering  action,  however,  some  degree  of 
signal  distortion  is  inevitable  unless  the  signal  is  a  constant  frequency  sinusoid.  Reducing  the  distortion 
to  a  tolerable  level  is  achieved  through  spectral  windowing;  i.e.,  to  appropriately  shape  and  flatten  (or 
level)  the  filter  response.  In  the  first  approach  attention  is  given  to  the  nature  and  characteristics  of  the 
signal  distortion  function.  The  trade-off  between  the  signal  and  SFT  parameters  is  defined  to  limit  the 
expected  degradation  to  a  tolerable  level.  In  the  revisited  approach  to  signal  transformation,  the  basic 
problem  of  signal  distortion  will  be  evaded  by  assuming  that  the  spectral  power  of  the  signal  is  bounded 
within  a  finite  section  of  the  signal  spectrum.  The  two  approaches  complement  one  another  and  yield 
an  insight  into  the  temporal  characteristics  of  the  transformation  algorithms. 

Signal  Description  and  Representation 

Over  an  extended  analysis  time  (approximately  T  seconds)  a  narrowband  signal  may  be 
represented  as 

«(/)  -  AU)e,<r0)  (la) 

where  the  phase  function  ♦  (f)  takes  the  form 

♦  (f)  -  2nlfei  +  /o  p(x)dx]  +  *0.  (lb) 

The  instantaneous  frequency  (or  inverse  wave-period)  of  the  signal  is  defined  as 

“-♦(f) -/,  +  *(f),  (1c) 
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where  fc  is  the  mean  frequency  and  v(t)  is  the  zero-mean  frequency  fluctuation  over  the  analysis  inter¬ 
val.  (A  "dot"  over  the  variable  is  used  to  denote  the  time  derivative.) 

In  the  analyses  to  follow,  two  subintervals  of  time  Tx  and  T2  are  to  be  employed  such  that  Tx  <  T2 
and  r,  is  much  much  less  than  the  extended  analysis  time. 

Letting  fcT2  *•  k0  +  (|8f|  <  1/2)  and  letting  r  ”  TJTX,  the  time  series  u(mTx)  becomes 

u(mT,)  -  A  -  eiUmke/,u0(mTx )  (2a) 

where 

u0(mTx)  -  A  («r1)/2’U/',*#<,"r,,+*01  (2b) 

and 

J  f 

g0(t)  —  +  ~  f0  v(x)dx.  (2c) 

The  sampled  time  series  u0(mTx)  represents  the  baseband  of  the  time  series  u(mT i)  referenced 
to  the  frequency  kJT2,  and  the  variable  g0(mTx)  is  a  running-time  average  (taken  at  mTx)  of  the 
instantaneous  frequency  deviation  (relative  to  baseband)  measured  in  units  of  1 IT2.  (The  merit  of  this 
form  of  notation  will  become  evident  when  we  consider  the  sectionalized  Fourier  Transform  of  the  sig¬ 
nal  u(r).) 

Sectionalized  Fourier  Transform 

Over  the  time  interval  T2  centered  at  mTx,  the  Sectionalized  Fourier  Transform  (SFT)  of  u(r)  is 
defined  as 

„  \  1  r'"rt+7VJ  — ,2*(*0+<.tf/r, 

Um(k0+n)-  -jrjmTrT2fi  uit)e  dt 

-  u0(«r,)e-fJ*""/'/_^j/fm(/)e'l,,|*"<n""lf/rjd//rJ  (3a) 

where 

/<„,(/)-  A(mTx  +  i)/A  (mTx)  (3b) 

and 

T  t 

*«(<)  “  *<•  +  ~  /„  v(mT\  +  x)dx.  (3c) 


(Care  must  be  exercised  in  the  interpretation  and  use  of  Eq.  (3b).  The  amplitude  function  Am(t )  is 
artificial  in  the  sense  that  A(mT\)  is  inserted  in  the  denominator  in  order  to  factor  out  u0(mTx)  in  Eq. 
(3a).  It  is  possible  therefore  for  ^(mf|)  to  be  zero,  in  which  case  Am(t)  makes  no  sense.  However, 
in  this  event  u0(mTx )  is  zero,  and  the  factor  Am(r)  should  rightfully  equal  only  the  numerator  term. 
For  many  practical  applications,  A(mTx  +  t)  will  be  essentially  constant  over  the  time  interval  T2, 
avoiding  the  possibility  of  singularities.) 

The  form  of  Eq.  (3a)  reveals  that  the  SFT  of  u(t)  yields  the  product  of  the  sampled  baseband  sig¬ 
nal  u0(mTt)  and  a  distortion  factor.  The  distortion  factor  is  a  function  of  the  spectral  selectivity  of  the 
SFT  and  the  static  and  dynamic  characteristics  of  the  signal.  Our  object  will  therefore  be  to  process  the 
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SFT  to  achieve  an  output  transform  which  approaches  u0(mT])  over  a  specified  signal  center-frequency 
and  bandwidth. 

Spectral  Windowing 

A  study  of  Eq.  (3)  suggests  the  use  of  a  spectral  window  comprised  of  J  sequential  frequency  bins 
n,  approximately  centered  at  the  frequency  k0/T2.  Therefore  let 

J-«  o"l 

Vn(k<£j)-  £  elll,mnl’Um(k0+  n) 

» — *0 


uQ(mT\)D„(gm\J), 


where 


Dm(gm-.J)  -  Am(t)e  -  2  e  *!Ti 


is  the  resulting  distortion  function  and  where  J  and  n0  are  chosen  to  essentially  constrain  the  signal 
power  within  the  spectral  window.  The  summation  within  the  above  integral  may  be  recognized  as  the 
Dirichlet  kernel  [21,22].  The  sum  reduces  to  [23] 

j-  i 

1+2  T  cos(2  irnt/Ti)  -  ^‘n  (when  J  is  odd) 

sm(ir//r2) 


±in  t/  T-y 

9  * 


]£cos  [ir(2n  -  Dr/rj  -  e 

n—  \ 


±ivt/T2  (w/r/T2) 

sin  (ni/Tj) 


(when  J  is  even). 


(The  exponential  factor,  when  J  is  even,  results  from  the  spectral  window  being  centered  midway 
between  two  spectral  bins  of  the  SFT.  The  sign  of  the  exponent  depends  on  whether  the  center  of  the 
window  is  located  one-half  bin-width  below  or  above  the  spectral  bin  k0 ) 


Using  the  Dirichlet  kernel  in  Eq.  (4b),  the  distortion  function  becomes 

«.<«/> - O.o> 


sin  (nt/Tj) 


dt/T,, 


provided  that  ±  1/2  is  added  to  the  parameter  5,.  (in  Eq.  (3c))  when  J  is  even.  (Assuming  that  v(t)  is 
symmetrically  distributed,  the  value  1/2  is  subtracted  when  bc  is  positive  and  added  when  fir  is  nega¬ 
tive.  This  procedure  is  required  if  the  spectral  window  is  to  most  efficiently  span  the  spectral  bandwidth 
of  the  signal.) 

Properties  of  the  Distortion  Function 

In  addition  to  the  window  parameter  J,  the  distortion  function  is  dependent  on  the  spectral 
characteristics  of  the  signal  u(t).  When  the  signal  dynamic  characteristics  are  sufficiently  slowly  vary¬ 
ing,  such  that  Ail)  and  v(i)  are  essentially  constant  over  time  intervals  of  r2  seconds,  the  distortion 
function  is  real  and  equal  to  the  spectral  window  function  Wj(xm).  That  is, 

-  Wj(xm)  -  Jo  cos  ( nxmt)dt .  (6a) 
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where  x*  ■■  8C  +  T-p(mT\)  when  /  is  odd  or  ic  T  0.5  +  Tp{mT\)  when  J  is  even.  By  expressing 
the  Dirichlet  kernel  by  its  equivalent  trigonometric  series  and  carrying  out  the  integration  prior  to  sum¬ 
ming,  Wj(x)  may  be  shown  to  be 


Wj(x) 


j- 1 

Y  sin  it  (x  -  n ) 
j-\  w(x-n) 
2~ 


sintrx 

irx 


(1 


j- 1 

+  2  X 


,2  - 


(6b) 


when  J  is  odd,  or 


Wjix) 


±(f-n 

L  sin  it  (x-  n ) 
,  w(x-n) 


1  2/2 

-  — cosjt(x±0.5)  X  (-D" 

IT  ^ 


n  -  0.5 


0-1 


(x  ±  0.5)2  -  (n  -  0.5)2 


(6c) 


when  /  is  even.  (The  shift  of  1/2  in  the  latter  relation  is  due  to  the  fact  that  the  center  of  the  spectral 
window  is  ±  one-half  bin-width  from  the  froth  bin.)  A  plot  of  Wj(x)  as  a  function  of  the  normalized 
frequency  is  displayed  in  Figs.  1  and  2  for  selected  odd  and  even  values  of  /,  respectively.  The  filter 
characteristics  are  somewhat  smoother  over  the  filter  passband  J/Tj  when  /  is  odd.  However,  as  /  gets 
larger,  the  difference  in  the  odd  vs  even  passband  characteristics  becomes  proportionately  smaller. 


Fi|  I  —  Passband  characteristics  of  the  spectral  window  function 
M^(jr)  for  /  odd 
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Fig.  2  —  Passband  characteristics  of  the  spectral  window  function 
Wt(x)  for  J  even 


In  addition  to  the  static  window  characteristics,  the  signal  distortion  is  also  a  function  of  the  signal 
dynamics.  The  expected  value  of  the  distortion  function  may  be  shown  to  approach  one  when  J  is  large 
and  the  spectral  bandwidth  of  the  signal  is  essentially  bound  within  the  filter  passband  (see  Appendix). 
Since  the  expected  value  of  the  distortion  function  was  computed  using  ensemble  averages,  it  will  prove 
informative  to  compute  this  function  over  the  time  interval  F2  for  a  representative  case.  Consider  then 
that  over  the  relevant  T2  time  interval,  the  amplitude  factor  AmU)  is  constant  and  the  frequency 
fluctuation  p(r)  can  be  represented  by  the  truncated  Taylor  series 

v(mT\  +  x)  •  v(mT\)  +  v(mT\)x.  (7a) 

From  Eq.  (3c)  and  the  definition  of  xm  (following  Eq.  (6a)) 

gmU)  -  xm  +  0.5Axm//r2  (7b) 

where 

Axm-p(mr,)r22.  (7c) 


(In  the  normalized  units  of  I/f2  Hz,  the  discrete  variable  xm  is  a  measure  of  the  instantaneous  fre¬ 
quency  deviation  from  the  window  center-frequency,  and  Axm  is  a  measure  of  the  change  in  instantane¬ 
ous  frequency  over  the  SFT  interval  T2.)  From  Eq.  (S)  then,  the  distortion  function  becomes 


sin  (it  Jill) 
sin  (nt/2) 


cos  {nxmt)dl. 


(8) 
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Except  when  Axm  is  zero  (in  which  case  the  distortion  function  reduces  to  the  spectra)  window 
function),  the  distortion  will  be  a  complex  number,  resulting  in  both  amplitude  and  phase  distortion. 
To  provide  an  indication  of  the  degree  of  distortion  that  can  exist,  Eq.  (8)  has  been  computed  as  a 
function  of  Ax„  (with  xm  -  0)  for  various  values  of  the  parameter  /  The  results  are  tabulated  in 
Table  1.  (The  ratio  hxJJ  is  the  fraction  of  the  spectral  window  over  which  the  instantaneous  fre¬ 
quency  varies  in  the  time  interval  T2.  Thus  when  xm  is  zero,  the  frequency  will  be  constrained  within 
the  spectral  window  for  bxm/J  less  than  one.)  The  data  in  Table  1  indicate  that  the  distortion  will  be 
relatively  minor  under  the  stipulated  conditions.  The  distortion  has  been  computed  for  values  of  xm 
other  than  zero  and  found  to  be  no  more  serious  than  that  shown  in  the  table,  as  long  as  the  instan¬ 
taneous  frequency  is  constrained  to  fall  within  the  bounds  of  the  spectral  window.  In  the  particular  case 
where  J  —  1  and  xm  —  0,  the  distortion  reduces  to 


A*<**;i)  -  (c(Vra7I)  +  «S(VU*J/2)1  (9) 

where  S(— )  and  C(— )  are  the  Fresnel  sine  and  cosine  integrals  of  the  indicated  argument  [24]. 


Table  1  —  Magnitude  and  Phase  of  Dm(gm\J )  as  a  Function  of 
Axm/7  for  the  Indicated  Values  of  J  (xm  —  0) 


J 

J  -  1 

3 

5 

7 

9 

0 

Magnitude 
Phase  (Deg.) 

1.000 

0.00 

1.000 

0.00 

1.000 

0.00 

1.000 

0.00 

1.000 

0.00 

±  0.2 

Magnitude 
Phase  (Deg.) 

0.999 

±3.00 

1.013 

4=1.81 

0.985 

±1.04 

1.016 

4=0.50 

0.986 

±0.16 

±  0.4 

Magnitude 
Phase  (Deg.) 

0.996 

±5.99 

1.051 

4=2.91 

0.949 

±0.49 

1.035 

±1.04 

0.988 

4=1.67 

±  0.6 

Magnitude 
Phase  (Deg.) 

0.990 

±8.98 

1.105 

4=2.85 

0.924 

4=2.56 

1.009 

±3.90 

1.044 

4=2.00 

±  0.8 

Magnitude 
Phase  (Deg.) 

0.983 

±11.96 

1.165 

=F1.53 

0.945 

4=7.17 

0.924 

±4.85 

1.085 

±2.42 

±  1.0 

Magnitude 

Phase  (Deg.) 

0.973 

±14.92 

1.224 

±0.94 

1.026 

4=10.71 

0.832 

4=0.25 

1.010 

±8.61 

SIGNAL  TRANSFORMATION  REVISITED 
Fourier  Series  Representation 

To  obuin  a  deeper  insight  into  the  sectionalized  Fourier  Transform  method  of  signal  filtering, 
basebanding,  and  sampling,  a  second  approach  will  be  taken  which  is  more  macroscopic  in  content. 
Utilizing  the  results  of  sampling  theory  [25,26],  it  has  been  shown  that  over  the  time  span  of 
mT\  -  Tj/2  <  t  mTx  +  T-Jl,  the  function  u(t)  may  be  expressed  by  the  Fourier  series 

u(t)  “  L  Um(k)ellwk,IT}.  (10) 

*  -  o 

where  Um(k)  is  the  sectionalized  Fourier  Transform  of  u(i)  (Eq.  (3))  taken  over  the  indicated  time 
interval.  Although  the  value  of  Um(k)  will  (in  general)  be  nonzero  over  all  k,  it  is  common  practice  to 
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set  bounds  on  the  range  of  k  over  which  )l/„(Jir)j  is  significant  for  certain  classes  of  signals.*  If  then 
u(t )  is  a  bandlimited  signal  whose  significant  spectral  energy  can  be  said  to  be  bounded  within  a  con¬ 
tiguous  sequence  of  J  spectral  bins  centered  at  approximately  k  -  fco.  a  suitable  approximation  of  uO) 
over  the  indicated  temporal  span  would  be 


*0  +  7-»0-t 

uO)  -  £  Um(k)eliwk,/Ti 

*  ~  *D_,,0 


or 


^i2wk0t/T} 


J-»  o-l 

I  tf„<*o  + 


"o 


n 


(lla) 


«„(/ )  -  e~'29k°,,T*u(t) 


-  Z  Um(ko+  n)et2’",/T\ 

«— «* 

Consequently,  from  Eqs.  (4a)  and  (lib), 

M^T-,)  -  K^Otoiy)  -  u0(mT\)Dm(g„\J). 


(lib) 

(12) 


The  above  relation  demonstrates  that  the  sectionalized  Fourier  Transform  can  serve  to  baseband, 
filter,  and  sample  a  bandlimited  function  without  basically  changing  the  temporal  characteristics  of  the 
function.  In  addition,  the  relative  simplicity  of  the  FFT  algorithm  permits  these  processing  operations 
to  be  performed  with  ease  on  modern  digital  computers.  The  factor  Dm(gm,J),  given  in  Eq.  (5),  pro¬ 
vides  a  measure  of  the  distortion  induced  by  the  process,  so  that  some  basis  is  available  for  the  selec¬ 
tion  of  the  parameters  7"2  and  J.  The  criterion  for  the  selection  of  the  sampling  rate  \/T\  is  also  well 
defined.  To  avoid  undersampling,  the  sampling  rate  should  be  equal  to  or  greater  than  the  (two-sided) 
filter  bandwidth  B}  —  J/Tj%  or  Tt  should  be  equal  to  or  less  than  the  Nyquist  interval  Tj/J  [27J. 


Equivalent  Temporal  Window 

The  temporal  counterpart  of  spectral  filtering  (or  windowing)  in  harmonic  analysis  is  to  shade  or 
to  weight  the  function  u(t)  over  the  time  window  mTt  -  Tj/2  <  t  <  mT{  +  T-J7  (28,291.  The  tem¬ 
poral  window  function,  equivalent  to  the  spectral  window  JP/Or),  is  simply  the  inverse  Fourier 
Transform  of  WjifT 2).  Therefore,  from  Eqs.  (6b)  and  (6c) 

wj(t)  -  r2J*~  lVj<JT2)ellwr'df 


7-1 


—  1  +  2  £cos(2w/»//7'2) 


sin  (irJt/Tj) 
sin(irf/r2) 


(I/I  <  772) 


(13a) 


‘Although  a  signal  cannot  be  both  band-  and  time-limited  in  a  pure  theoretical  sense,  this  concept  has  proven  quite  useful  in 
practical  applications.  An  excellent  discussion  of  the  problem  is  found  in  Ref.  21.  pp.  121-132. 
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for  J  odd,  and 


wj(t)  -  e±in',T 2  2  £cosbr(2/i  -  Dz/T}] 


*nrt/r,  sin  (irJt/T2) 
sin(ir//r2) 


(I/I  <  TJ 2). 


for  J  even.  (The  complex  exponential  factor  for  J  even  is  a  consequence  of  the  spectral  window  being 
centered  midway  between  two  spectral  bins.)  A  plot  of  the  temporal  window  function  for  various 
values  of  J  is  shown  in  Fig.  3.  For  J  —  1,  the  weighting  is  constant.  For  7-2,  the  weighting  magni¬ 
tude  is  a  simple  cosine  function.  As  7  becomes  larger,  the  temporal  window  (Dirichlet  kernel)  more 
closely  approaches  a  (sin  x)(x  function  over  the  interval  T2. 


tm(r  Jt/T2) 

*J{,>  * 


7  —  Bjltl 


Fig.  3  —  Temporal  window  function  wj{i)  for  selected  values  of  J 


If  now,  over  the  time  interval  mf,  -  TJ 2  ^  f  <  mf,  +  TJ2,  one  defines  a  new  function  v(r), 

v(/)  -  wAt  -  mT,)u(t),  (U) 


W' 
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its  sectionaiized  Fourier  Transform  can  readily  be  shown  to  be 


1  p»T,+rj/2 

f2  J«r,-V2 


wj(t  —  mT\)u(t)e  '2"k,ITi(/t 


where 


J-n  ,,-i 

-  £  elUmn/rUm(k  +  n)  -  Km(*;/), 

"—"o 


IU*) 


_L  f 

r, 


mTl  +  T}/2 
■772 


u(t)e 


— itorktlT-} 


it 


(15a) 


(15b) 


and  where  n0  is  (/—  l)/2  when  /is  odd,  and  //2  or  //2  —  1  when  /  is  even  (depending  on  the  sign  of 
the  exponent  in  Eq.  (13b)).  Equation  (15a)  is  identical  to  Eq.  (4a),  except  that  here  the  spectral  bin  k 
is  general  and  need  not  be  restricted  to  fitting  any  particular  center  frequency.  Thus  the  spectral  win¬ 
dow  may  be  incremented  across  the  frequency  band  in  steps  of  one  to  /  bin-widths.  Computationally,  it 
is  more  economical  to  implement  the  window  function  in  the  spectral  domain  when  only  one  or  a  few 
spectral  bands  need  be  examined  for  signals.  There  is  also  merit  in  the  fact  that  k  may  be  indexed  in 
less  than  /  increments  to  smooth  or  effectively  eliminate  any  'picket-fence*  or  "scallop"  effect  between 
windows  [28]. 

A  study  of  the  temporal  window  function  (Fig.  3)  reveals  that  the  significance  of  u(r),  in  the  for¬ 
mation  of  its  Fourier  Transform,  decreases  rapidly  as  /  deviates  from  ±  T2/2J.  That  is,  as  |/|  becomes 
greater  than  1/2 Bj,  the  weight  given  to  u(/)  becomes  appreciably  reduced,  so  that  its  significance  in  the 
construction  of  the  resulting  Fourier  Transform  Um(k)  is  reduced.  This  is  why  the  value  of  Axm,  in 
the  distortion  function  Dm(gm\J),  can  become  proportionally  larger  with  /  without  seriously  altering  the 
transform  characteristics  (see  Table  1).  The  restrictions  on  the  rate  of  change  in  the  amplitude  function 
A  (r)  is  also  proportionally  reduced.  Another  way  of  looking  at  it  is  that  as  /  becomes  larger,  the 
bandwidth  of  the  spectral  window  increases  (assuming  T2  remains  fixed).  And  consequently  the  signal 
dynamics  can  be  correspondingly  more  rapid  (spreading  the  power  spectral  density  of  the  signal)  to  fill 
the  wider  window,  without  seriously  degrading  the  transformed  output. 


Doppler-Induced  Distortion 


When  a  signal  source  is  in  motion  in  a  transmission  medium,  the  spectral  energy  of  the  signal 
suffers  a  Doppler  shift.  The  effect  of  the  Doppler  shift  is  to  compress  (or  expand)  the  time  scale  of  the 
original  signal  [30].  Thus,  the  signal  u(t)  is  transformed  into  u{(l  +a0)/}  where  a0  (a0  «  1)  is 
known  as  the  time  scale-factor  shift  or  Doppler  ratio.  From  Eqs.  (1)  and  (2)  it  is  easy  to  show  that 


«0{(1  +  a0)m7’||  -  A2mu0(mTl)e'2,rma'>lrc+‘’i,mT'))r\  (16a) 

where  u0[(l  +  and  u0(mT,)  are  the  respective  band-shifted  signals  (relative  to  ko/T2  Hz),  and 

where 


and 


A2m  -  A {(1+  a0)mT]}/A (mTt)  (16b) 

v2(mT\)  =  v{(l  +  ao/2)mT\).  (16c) 


When  a0  is  sufficiently  small  so  that  rm0  <  r/J,  A2  m  will  be  close  to  unity.  (The  amplitude  distortion 
factor  A2m  is  due  to  the  time  compression  of  the  Doppler-shifted  signal.)  For  the  purposes  of  this 
paper,  A  (/)  is  considered  to  vary  sufficiently  slowly  so  that  insignificant  error  will  result  in  assuming 
that  /d2.m  —  1. 
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Since  the  effect  of  the  time  compression  is  to  slightly  shift  the  spectral  power  of  the  signal  u(t) 
into  a  new  band,  it  may  be  desirable  to  translate  the  Fourier  Transform  sequence  (spectral  window)  to 
accommodate  this  Doppler  shift.  To  optimumly  accomplish  this,  let 

aofcTt  -  na  +  «,  (where  -1/2  <  c  <  1/2).  (17a) 

Following  the  procedure  given  in  Eq.  (11),  it  may  be  verified  that 

+«o)«7'i)  -  Vi ,*(*o  +  na\J) 

J-i0- 1 

-  I  eiUm”''U2,m(k0+  na  +  n) 


~  fi0(/n7'1)e,J*,",,/,+aor|',J<"r|)l,  (17b) 

where  U2m(k)  is  the  sectionalized  Fourier  Transform  of  u{(l  +  a0)r).  The  translation  of  the  Fourier 
Transform  sequence  by  n„  has  the  effect  of  centering  the  spectral  window  on  the  Doppler-translated 
signal  spectrum  to  minimize  the  output  signal  distortion.  (If  the  spectral  window  is  sufficiently  broad  to 
adequately  encompass  the  Doppler  shift,  this  step  would  be  unnecessary.)  Another  way  of  looking  at  it 
is  that  the  far  left-hand  side  of  Eq.  (17b)  is  the  baseband  for  the  Doppler-shifted  signal.  Although  a 
frequency  translation  of  k<J  T2  represents  baseband  for  the  signal  u  (r),  a  frequency  translation  of 
(A:0  +  na)/T2  is  required  to  represent  baseband  for  the  Doppler-shifted  signal  k((1  +  a0)r). 

For  our  purposes  we  shall  assume  that  the  bandwidth  of  the  spectral  window  is  sufficient  to  ignore 
the  Fourier  Transform  distortion  factor  (permitting  us  to  drop  the  'tilde"  from  the  functional  relations). 
And  we  shall  assume  that  ma0  is  sufficiently  small  to  ignore  the  amplitude  distortion  factor  A2„.  (We 
shall  later  develop  the  restrictions  on  a0  and  the  analysis  time  to  permit  this  realization.)  The 
significant  relations  relative  to  Doppler-shifted  signals  are  then 

Mold  +  a0)«r,)  -  V2.m{ko,J)  -  V2.Jk0  +  na,J) 

_e‘»m*'ScT,yjko (lg) 

The  first  exponential  factor  in  the  right-hand  side  of  Eq.  (18)  reveals  that  a  Doppler  shift  produces  a 

linearly  varying  phase  rotation  on  the  original  signal.  The  rate  of  phase  rotation  is  proportional  to  the 

product  of  <*0  and  the  mean  signal  frequency  fc.  (For  a  cw  signal  this  will  be  the  only  phase  distortion.) 
However,  the  zero-mean  fluctuating  frequency  v{t)  introduces  a  nonlinear  phase-shift  which  must  also 
be  taken  into  consideration.  The  degradation  effect  of  this  latter  factor  will  be  addressed  in  a  latter  sec¬ 
tion  of  the  paper. 

Effect  of  Time  Shifts 

Consider  now  the  effect  of  a  simple  time  translation  r0  on  the  signal  u(t).  Letting  r0  — 
(mo+«i)T|  where  -1/2  <  e,  <  1/2,  it  is  easy  to  show  that 

u0(mT,  +  r0)  -  ^i,mM0{(«  +  m0)7',}e'2',,'^+,'',('B+"4,r,,,r'.  (19a) 

where  u0(«r1  +  r0)  and  u0{(m  +  m0)7'|}  are  the  respective  basebanded  signals,  and  where 

A  ii(W  -  A  (mT|  +  r0)/A  {(m  +  mo)  T|J  =*  1  (19b) 

and 

v,(mr,)-vl(m  +«,/2)r,).  (19c) 
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(Since  the  amplitude  (unctions  in  Eq.  (19b)  differ  in  time  by  less  than  1/2 Tu  negligible  error  will  result 
in  assuming  that  Alm  -  1.) 

Again,  following  the  procedure  given  in  Eq.  (11),  it  may  be  verified  that 

y-*o-i 

fl,(«r,  +  To>-  K, ..(**/>-  I  e^'U^Jko+n) 

«—  "0 

-  ^2™,l/r+,1l<m+m0>rIll (20) 

where  U\  m(k)  is  the  sectionalized  Fourier  Transform  of  u(t  +  t0). 

Assuming  that  the  bandwidth  of  the  spectral  window  is  sufficient  to  ignore  the  Fourier  Transform 
distortion  factor  (Eq.  (12)),  the  "tilde"  may  be  dropped  from  the  above  relation.  Consequently,  the 
effect  of  the  residual  parameter  is  to  cause  a  fixed  phase-shift  (2vt\fcT\)  and  a  fluctuating  phase- 
shift  (due  to  the  fluctuating  frequency  v(r))  on  a  signal  that  would  be  delayed  an  even  multiple  of  T\ 
seconds. 

DESIGN  CONSIDERATIONS 
Noise  Power  Output 

Since  the  signal  channel  will  generally  be  contaminated  by  broadband  noise,  it  will  be  of  interest 
to  determine  the  noise  output  of  the  filter  W(JT-^  inherent  in  the  sectionalized  Fourier  Transform. 
Assuming  a  broadband  noise  power  spectral  density  of  N0  watts  per  Hz,  the  accumulated  noise  power 
over  the  spectral  window  is  (employing  Parseval’s  theorem;  see  Ref.  25,  p.  65) 


PN- N0  J”  W}(JT2)df  -  ^r2- J#"  W}(x)dx 


r  sin2(ir/f)  . 
T2  **°  sin2  (it/) 


-  No  JlT-i  -  N0Bj, 


(21) 


where  Bj  -  J/T2  is  the  bandwidth  of  the  spectral  window  (see  Figs.  I  and  2).  (It  is  also  the  equivalent 
noise  bandwidth  of  the  window  (281.) 


Output  Signal-to-Noise  Ratio 


If  py{f)  is  the  signal  power  spectral  density  function  and  /  is  frequency,  measured  relative  to  the 
center  of  the  spectral  window,  the  output  signal  power  will  be 

Pu  -  J_^(/>»7(/T2)«//.  (22) 


When  the  signal  power  is  uniformly  distributed  over  the  band  f\  to  f2,  the  output  signal  power 
becomes 


Pu 


<uHt)> 


Zjiyj)  -  Zj(y i) 


(23a) 
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where  y  —  2 f/Bj  and  where 

Zjiy)  -  WjUx/Ddx. 


(23b) 


The  function  Zjiy)  is  tabulated  in  Table  2  for  values  of  y  lying  within  the  main  lobe  of  the  spectral 
window.  From  Eqs.  (21)  and  (23)  then,  the  output  signal-to-noise  power  ratio  becomes 


PJPn 


<u1(t)>  Zj(y 2)  ~  Zj(y  1) 
N0Bj  y2  -  y. 


(24a) 


And,  if  the  signal  power  spectral  density  is  symmetrically  located  within  the  filter  window,  the  relation 
reduces  to 


PJPn 


<u2(t)>  Zj(y) 
NoBj  y 


(24b) 


For  a  properly  designed  system,  Bj  should  be  chosen  to  efficiently  contain  the  signal  power 
without  significant  excess  or  spillover.  Ideally,  y  should  be  one  to  avoid  an  excess  of  unwanted  noise 
and  interference.  From  Table  2  (and  y  —  1),  it  is  seen  that  for  J  greater  than  one  the  output  signal- 
to-noise  ratio  is  very  nearly  ideal.  Even  for  J  equal  to  one,  the  loss  in  sign 'll -to-noise  ratio  is  only 
slightly  in  excess  of  one  decibel.  Consequently,  from  the  standpoint  of  signal-to-noise  ratio,  there  is  no 
strong  motivation  for  choosing  a  value  of  J  greater  than  one. 


Table  2  —  Integral  of  the  Square  of  the  Spectral  Window  Function 
WjUxI 2)  over  the  Limits  0  to  y  (see  Eq.  (23b)) 


1 

2 

3 

■ 

5 

6 

■■ 

8 

9 

_ 

10 

0.0997 

0.1611 

kin  dill 

kjXjMJJ 

0.1186 

0.1006 

kz;  1 

ijijua 

mffl 

0.1491 

0.1727 

0.1516 

0.1978 

0.1567 

0.2219 

ffifn 

OL  j 

0.2458 

0.3900 

0.2043 

0.2664 

0.2513 

0.2444 

0.30 

0.2927 

0.4602 

0.3123 

0.2909 

0.3046 

0.3095 

0.2994 

0.2919 

0.35 

0.3385 

0.5266 

0.3387 

0.3468 

0.3533 

0.3646 

0.3518 

0.3446 

0.40 

Lt 

0.5886 

0.3774 

0.3866 

0.4003 

0.4146 

0.4042 

0.4014 

0.45 

0.6459 

EGSI 

0.4443 

0.4461 

0.4593 

0.4551 

0.4585 

0.50 

0.4674 

0.6984 

0.5409 

0.5140 

0.4761 

0.4918 

0.5005 

0.5034 

0.5115 

0.55 

0.5070 

0.7459 

0.5978 

mm 

ji 

0.5294 

0.5411 

0.5492 

0.5584 

0.60 

0.5448 

0.7884 

0.6531 

0.6546 

0.5893 

0.5941 

SMS*!* 

0.65 

0.5807 

0.8259 

0.7059 

E2 

yglll 

0.6440 

0.6406 

0.6419 

0.70 

0.6145 

0.8586 

0.7552 

0.7300 

0.7229 

0.7030 

0.6931 

0.6916 

0.6879 

0.75 

0.6463 

0.8867 

0.8003 

0.8361 

0.7881 

0.7901 

0.7649 

0.7581 

PETTI 

0.7428 

0.80 

0.6760 

0.9106 

0.8403 

0.8414 

0.8516 

0.8261 

0.8065 

0.8S 

0.7036 

0.9305 

0.8749 

0.9187 

0.8874 

0.9034 

0.8820 

0.8874 

■  ,*\9 

0.8729 

0.90 

0.7291 

0.9467 

0.9038 

0.9470 

0.9245 

0.9430 

0.9278 

0.9373 

0.9257 

0.9307 

0.95 

0.7524 

0.9598 

0.9269 

0.9673 

0.9520 

0.9700 

0.9605 

0.9640 

0.9707 

1.00 

0.7737 

0.9701 

0.9447 

0.9808 

0.9703 

0.9859 

0.9801 

0.9888 

0.9851 

0.9908 

1  +  1/J 

0.9028 

0.9949 

0.9762 

0.9959 

0.9876 

0.9967 

0.9919 

0.9972 

0.9941 

0.9976 

Spectral  Leakage 

An  important  facet  in  the  design  of  the  spectral  window  (selection  of  the  parameter  J  )  is  the 
spectral  leakage  resulting  from  the  window  sidelobes  (28J.  These  sidelobes  will  cause  signals  remote 
from  the  window  bandwidth  to  appear  at  the  filter  output,  even  though  they  are  attenuated.  And  if  the 
remote  signals  are  sufficiently  strong,  they  can  seriously  interfere  with  those  signals  falling  within  the 
filter  window.  Consequently,  it  is  important  that  the  filter  window  reject  those  signals  whose  spectral 
energy  falls  outside  of  the  spectral  window  to  the  degree  that  is  practical. 

From  Eq.  (6),  the  magnitude  of  the  spectral  window  sidelobes  is  closely  approximated  by 


Wj(  —  -t-  m)  | 


(/  +  2m)2 


when  /  is  odd,  or 


1  w-{LT-  +  ”11  ■  71 1  (-'),o7I#rj£z-iSi| 

when  J  is  even,  where  m  —  1,2,  ...  is  the  sidelobe  index  along  the  frequency  axis. 


For  large  values  of  m  (J  «  m ), 
I  WAj  +  m)  | 


VOTu*  -  TTr  <f°' J  M 


W0(~~y~L  +  m)  | 


ir(/  +  2m)2  lirJf2 


(for  J  even). 


where  /  is  frequency  measured  relative  to  the  center  of  the  spectral  window.  For  Bj  constant,  it  is  seen 
that  the  magnitude  of  the  remote  filter  sidelobes  is  inversely  proportional  to  J.  A  more  interesting  fact 
is  that  the  magnitude  of  the  sidelobes  decays  at  a  rate  of  12  dB  per  octave  of  frequency  when  /  is  even, 
and  only  6  dB  per  octave  of  frequency  when  J  is  odd.  Thus,  from  the  standpoint  of  interference  rejec¬ 
tion  it  will  be  more  productive  to  make  /  even. 

Plots  of  the  spectral  window  characteristics  for  /  odd  and  even  are  shown  in  Figs.  4  and  S.  The 
frequency  axis  of  the  curves  is  scaled  in  units  of  the  window  bandwidth  in  each  case  for  comparison 
purposes.  The  advantage  of  making  /  even  in  lieu  of  odd  is  quite  apparent.  Further,  since  the  side- 
lobe  density  is  J  lobes  per  window  bandwidth,  interferring  signals  whose  spectral  power  is  spread  over 
one  or  more  sidelobes  will  be  attenuated  approximately  4  dB  below  the  indicated  sidelobe  envelope. 

CORRELATION  PROCESSING 

Magnitude-Squared  Coherence 

The  magnitude-squared  coherence  function  (MSC)  of  two  signals  st(/)  and  s2(t)  is  defined  as 

rV) _ IgiML  (27) 

J  St, (/>$«(/> ' 


where  Sn(f)  and  S2i(f)  are  the  power  spectral  densities  of  S|(/)  and  sj(f)  respectively,  and  S)2(/)  is 
the  cross-power  spectral  density  (Fourier  Transform  of  the  cross-correlation  function).  Note  that  the 
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Fig.  4  —  Sidelobe  leakage  characteristics  of  the  spectral  window  function  IF/I.e)  for  odd  values  of  the  parameter  J 
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FREQUENCY  IN  UNITS  OF  WINDOW  BANDWIDTH  J/Tj 


Fig.  5  —  Sidelobc  leakage  characteristics  of  ihc  spcciral  window  function  Wylvl  for  even  values  of  the  parameter  J 
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MSC  will  range  between  zero  and  one  depending  upon  the  magnitude  of  the  cross-power  spectral  den¬ 
sity. 

Perhaps  the  most  unique  property  of  the  MSC  is  its  invariance  under  linear  operations.  That  is,  if 
V|(r)  and  v2(/)  are  the  result  of  linear  time-invariant  operations  on  jt(/)  and  s2(t),  then  the  MSC  of 
V|(r)  and  v3(/)  will  be  identical  to  the  MSC  of  S|(r)  and  s2(t)  {4}.  This  will  be  true  even  though  the 
correlation  coefficient  between  v, (r)  and  v20)  may  differ  radically  from  the  correlation  coefficient 
between  s\(t)  and  s2(t).  Thus,  the  square  root  of  the  MSC  and  the  correlation  coefficient,  although 
somewhat  related,  are  truly  different  concepts. 

MSC  Estimate 


The  magnitude-squared  coherence  estimate  has  been  defined  as  [4] 


T2U) 


I  <  S{m(k)S2'jk)  >  p 


2,m  ' 


<  |S,.„(/c)P  >  <  \S2Ak)\2  > 


(28) 


where  the  indicated  average  is  computed  over  a  given  analysis  interval,  and  Sm(k)  is  the  sectionalized 
Fourier  Transform  of  the  relevant  temporal  function  computed  over  time  T2.  (The  asterisk  denotes  the 
complex  conjugate.)  In  the  notation  above,  the  time  index  is  m  (t  ■  mT\^ mT2)  and  the  frequency 
index  is  k  (/*  -  k/T2).  The  numerator  in  the  above  equation  forms  the  estimate  of  the  cross-power 
spectral  density  between  the  two  relevant  signals  spr)  and  s2(t). 

The  MSC  estimate  has  received  considerable  attention  in  the  literature  as  a  sample  test  statistic 
for  coherence  estimates  14-16].  However,  depending  upon  the  length  of  the  analysis  interval  and 
several  other  factors,  the  estimate  may  not  be  a  good  estimate  of  the  MSC.  One  can  readily  perceive 
that  the  estimate  is  no  longer  invariant  with  linear  operations  on  the  two  signals,  but  can  vary  appreci¬ 
ably  depending  on  the  nature  of  these  operations.  In  fact  (as  shall  be  subsequently  shown),  when  the 
spectral  power  of  the  two  temporal  signals  is  essentially  bounded  within  the  spectral  bin-width  of  the 
Fourier  Transform,  the  square  root  of  the  MSC  estimate  closely  approximates  the  normalized  envelope 
correlation  function  of  the  two  signals. 


Normalized  Correlation  Envelope 

Consider  now  that  S|(f)  and  s2(r)  are  two  real  narrowband  signals  present  at  two  sensors.  The 
normalized  two-dimensional  correlation  function  (NC)  of  the  signals  (over  an  extended  analysis  inter¬ 
val)  is  defined  as  131) 


■y  (r ,  a ) 


<  S\(t  —  r)sj(r  +  at)  > 

V  <  5\  (r  -  r)>  < s2  (/  +  at)  > 


(29) 


where  the  indicated  averaging  is  carried  out  over  the  analysis  interval.  The  resulting  NC  can  generally 
be  written  as  the  product  of  a  slowly  varying  correlation  envelope  (NCE)  function  *(T,a),  and  a 
sinusoidal  carrier  function  C(r,a)  [31).  By  repeating  the  NC  with  r  shifted  by  one-quarter  of  a  cycle 
of  the  carrier  frequency,  the  resulting  NC  will  be  in  quadrature  with  the  original  NC.  (The  minute  shift 
in  r  will  not  significantly  change  the  value  of  the  correlation  envelope  function.)  Thus,  the  NCE  may 
readily  be  computed  as  the  square  root  of  the  sum  of  the  squares  of  the  NC  and  the  quadrature  NC. 

The  indicated  procedure  for  determining  the  NCE  is  computationally  awkward  and  inefficient. 
Further,  since  the  desired  signals  are  generally  contaminated  by  noise  and  interfering  signals,  some 
form  of  filtering  is  desired  around  the  relevant  signals  to  improve  their  signal-to-background  ratio  prior 
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to  correlation.  The  technique  of  basebanding  and  low-pass-filtering  the  relevant  signals  prior  to  correla¬ 
tion  is  a  suitable  alternative,  but  this  is  computationally  unattractive  if  carried  out  in  the  normal  step¬ 
wise  fashion.  However,  the  earlier  analyses  have  shown  that  basebanding  and  spectral  filtering  may  be 
accomplished  quite  readily  through  the  use  of  the  sectionalized  Fourier  Transform.  This  suggests  that 
the  following  algorithm  may  serve  as  a  convenient  estimate  of  the  NCE. 


\<em»'vx.m-m(kJ)vUkM>\ 
/< I < I v2.m(k-,J) IV 


where  the  symbology 


i  M/2-1 

<(— )>  implies  —  £  (• 

M  m— M/2 


and  where  mr  (an  integer)  reflects  the  time  displacement  (tau  variable),  and  q  (an  integer)  reflects  the 
time  scale-factor  shift  or  Doppler  ratio  (alpha  variable)  between  the  two  signals.  The  transforms  in  the 
relation  are  defined  by  Eq.  (IS),  with  subscripts  added  to  denote  the  signals  within  the  two  Ar-bin  chan¬ 
nels  being  processed.  The  exponential  factor  in  the  numerator  of  the  above  relation  serves  as  the 
ambiguity  kernel  exp(/2ira/r)  [32]  to  Doppler-shift  the  transform  V2m(k\J).  This  is  suggested  from 
Eq.  (18)  as  a  method  of  (approximately)  compensating  for  any  time  scale-factor  shift  between  the  sig¬ 
nals  in  the  two  channels,  t 

The  optimum  choice  for  T\  in  the  NCE  estimate  is  the  Nyquist  interval  T-JJ  —  1/fly 
( r  -  TJT\  “  7).  With  this  choice,  the  total  analysis  time  is  A/F,  and  the  correlation  integration  time 
is 

F-  (A/-  1)F,  -  (A/-l)/fl;  (31a) 

or 

A/-flyF+l  (31b) 

An  explicit  expression  for  the  r  variable  is,  of  course, 

r  -  mrT\  —  mj Bj  ;  (32) 

however,  the  explicit  Doppler  ratio  for  a  given  signal  will  depend  on  the  mean  frequency  of  that  signal 
over  the  processor  analysis  interval  (Eq.  (16a)).  A  suitable  estimate  of  the  Doppler  ratio  is  given  by 

*  '  2M*  *  2*T  <33,) 

where  /*  “  k/T2  is  the  approximate  center-frequency  of  the  spectral  window.  When  k  is  optimum  for 
a  given  signal  (viz.,  when  the  spectral  window  is  most  nearly  centered  about  the  signal  spectrum),  the 
error  in  the  estimate  will  be 

I— I  -  «  rrr  «»> 

a  k  27/* 

where  ftr  (|6C  |  <  1/2)  is  the  difference  between  fc  and  the  nearest  harmonic  k/T2.  The  choice  of  the 
index  q  to  vary  the  Doppler  ratio  in  increments  of  approximately  l/2/rF  was  made  to  limit  the  ’picket- 
fence*  or  ’scallop”  loss  in  the  a-dimension  to  less  than  1  dB  [281. 

tR.D.  True  blood  of  the  Navel  Ocean  Systems  Center  (NOSC)  employed  this  technique  in  connection  with  the  magnitude- 
squared  coherence  (MSC)  estimate  in  the  early  I970's.  The  results  of  his  (and  subsequently  other's)  investigations  indicated  that 
the  compensation  was  adequate  for  the  signals  and  parameters  used  in  the  MSC  estimate. 
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Computation  of  NCE  Estimate 

Computing  the  two-dimensional  NCE  estimate  must  certainly  be  considered  a  formidable  task  as  a 
result  of  the  variables  mT  and  q  over  which  a  correlation  surface  is  mapped.  In  practice,  the  number  of 
points  required  to  satisfactorily  map  the  surface  can  range  from  several  thousand  to  the  tens  of 
thousands.  However,  the  present  state-of-the-art  in  computer  technology  (using  modern  array  proces¬ 
sors)  is  such  that  adequate  computational  speed  is  available,  providing  the  processing  algorithms  are 
suitably  matched  to  the  computer  architecture.  In  the  case  of  the  transforms  Km  (*;/),  these  are  ideally 
suited  for  array  processing  using  conventional  DFT  (discrete  Fourier  Transform)  formats.  One  may 
also  perceive  that  the  numerator  of  the  NCE  estimate  (Eq.  (30))  has  the  form  of  a  DFT  relative  to  the 
Doppler-ratio  variable  q.  This  considerably  reduces  the  required  programming  and  the  computational 
time  in  the  realization  of  a  NCE  surface. 


The  variation  in  the  r-dimension,  although  not  as  expedient  to  reproduce,  involves  incremental 
translations  of  the  one  transform  Vm(k\J)  relative  to  the  other.  When  the  equivalent  bandwidth  of  the 
desired  signal  spectral  power  is  greater  than  one-half  the  width  of  the  spectral  window,  the  scallop  loss 
(due  to  discrete  sampling  in  the  r-dimension)  can  exceed  1  dB  and  may  approach  4  dB  (when  the 
equivalent  bandwidth  is  equal  to  the  width  of  the  spectral  window).  Under  these  circumstances,  the 
scallop  loss  can  be  reduced  by  computing  a  second  series  of  Vm(k\J)  (for  one  of  the  signals),  which  is 
temporally  interlaced  with  the  original  set.  This  is  equivalent  to  choosing  Tt  to  be  1/2 Bj  for  the  one 
signal,  and  then  sequentially  interlacing  the  odd  and  even  sets  of  resulting  sectionalized  Fourier 
Transforms  in  computing  the  NCE  estimate. 


Although  somewhat  complicated,  the  indicated  techniques  (or  modified  versions  thereof)  for 
computing  rather  extensive  ambiguity  surfaces  have  been  accomplished  with  relative  ease  on  suitable 
array  processors  within  the  past  decade  (18-20,311. 

Validation  of  the  NCE  Estimate 


The  algorithm  for  estimating  the  NCE  (given  by  Eq.  (30))  can  be  applied  in  two  ways.  First,  it 
can  be  used  to  study  the  ambiguity  surface  features  for  a  broad  class  of  temporal  functions.  And 
second,  it  can  be  used  to  detect  (and  estimate  the  parameters  for)  common  signals  which  differ  in  time 
alignment  and/or  Doppler  ratio. 

In  the  first  application,  Vl  m(k\J)  —  Vi  m(k\J)  and  the  parameters  q  and  m,  are  used  to  map  the 
autocorrelation  envelope  over  the  ra -plane.  Letting  the  spectral  window  be  centered  on  (and  encom¬ 
pass)  the  spectral  energy  of  the  signal  u(r),  it  is  seen  that  Eq.  (30)  closely  approximates  (from  Eqs. 
(18)  and  (20)) 


*  JT'a)  " 


|<«0(mri  -  r)i/o*Kl  -a)mT\)e 

y/ <  -  t)(^>  <  |m0{0  +  a)mTi\i> 


where 


t  —  mr/Bj  and  a  —  —qBj/2Mfc. 


(34a) 


(34b) 


When  both  mT  and  q  are  zero,  the  NCE  estimate  equals  unity.  Further,  along  the  r  axis  (a  —  0),  the 
estimate  of  the  NCE  is  essentially  precise.  However,  when  a  is  nonzero,  the  estimate  is  degraded  as  a 
consequence  of  the  nonlinear  phase-distortion  factor.  (The  degree  that  this  factor  influences  the  NCE 
estimate  will  be  determined  shortly.) 

In  the  second  application,  consider  that  the  common  signal  is  uU)  (see  Eq.  (1))  and  that  st(r)  — 
w(r  4-  r0)  and  s2(r)  •  u(t  +  a0t).  In  this  application,  the  parameters  r0  and  a0  are  unknown  and  will 
need  to  be  estimated  from  the  NCE  function  topology.  (The  NCE  estimate  is  expected  to  peak  at  the 
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point  where  r  and  a  compensate  for  the  parameters  r0  and  a0  respectively.)  In  practice,  mr  and  q  are 
systematically  sequenced  over  a  set  of  values  which  will  encompass  the  anticipated  range  of  r0  and  a0. 
Since  the  intent  of  the  processing  is  to  compensate  for  time  and  Doppler  differences  in  the  received  sig¬ 
nals,  it  will  be  advantageous  to  shift  the  spectral  window  periodically  as  discussed  in  relation  to  Eq. 
(17).  (The  shift  in  the  spectral  window  function  is  to  ensure  that  the  two  spectral  windows  span 
approximately  the  same  portion  of  the  signal  spectrum  in  the  common  signal.) 

An  algorithm  for  determining  the  bin-shift  parameter  is  developed  as  follows.  With  M  (see  Eq. 
(30b))  chosen  so  that  M/J  is  an  integer,  let  q'  be  a  modulo  integer  defined  as 

q'  -  (q  -  JM/J)  Mod  (2 jM/J)  -  jM/J  (35a) 

where  J  (j  <  J)  is  the  desired  number  of  Fourier  Transform  bins  to  be  translated  with  each  modulo 
sequence.  (When  J  is  small,  it  may  be  expedient  to  make  J  ™  1;  however,  when  J  is  large,  computa¬ 
tion  time  can  be  saved  by  translating  the  spectral  window  in  larger  increments  without  seriously  degrad¬ 
ing  the  sensitivity  of  the  correlation  processor.)  The  shift  parameter  na  is  then 

na  -  (q  -  q')J/2M,  (35b) 

or,  in  terms  of  na,  the  value  of  q  is 

q  -  —  na  +  q.  (35c) 


From  Eqs.  (12),  (18),  and  (20)  then,  it  may  be  verified  (with  a  little  algebraic  manipulation)  that 
**«  T’“  y/<  I  \*>  <  |  VJm(k0  +  naxJ)  |*> 


l<u0(mr,  -  t')mq*((1  +  a')mTl)e~nwmia,f}imr')  ~  | 

>/< -  t')I2><|m0I(1  +  o')m7*,}|2> 

where 

r’  -  r  -  r0  -  mJBj  -  t0 
a'  “  <*o  -  a  «  a0  —  qBj/2Mfc 
vjlmTf)  -  i»{(l  4-  a<J2)mTi) 

and 

v'i(mT\)  -  vj(l  +  a'/2)mr|). 


(36a) 


(36b) 

(36c) 

(36d) 

(36e) 


The  above  relation  verifies  that  the  NCE  estimate  is  an  accurate  representation  of  the  NCE, 
except  for  the  nonlinear  phase-distortion  factor  in  the  numerator.  Of  course,  when  a  -  0  (q  —  0  and 
a'  —  a0)  the  estimate  is  precise,  provided  the  spectral  energy  of  both  signals  falls  within  the  spectral 
window.  (This  results  from  the  fact  that  no  Doppler  compensation  has  been  employed.  The  nonlinear 
phase-distortion  is  due  solely  to  the  convenient  method  chosen  for  Doppler  compensation.)  When 
o  is  equal  to  (-1/2  <  «,  <  1/2)  and  a0-qB,/2Mfc  is  equal  to  t:Bj/2Mfc 

(-1/2  <  <}  <  1/2),  the  estimate  will  maximize  and  equals 


X.  (r.or) 
*0 


JW«| 

I  <Rme 


(37a) 
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where 

Rm  -  \u0(mTl)\2/<\u0(mTi)\2>  -  A2(mT{)/ < A2(mT,)>  (37b) 

r  -  t0  -  (JBj  (-1/2  <  «,  <  1/2)  (37c) 

a  -  «0  “  t:Bj/2Mfc  (“1/2  <  «2  <  1/2)  (37d) 

i-,(mr1)-^(ffl+t,/2)r,|  (37e) 

and 

vi(mT\)  -  */{(!  +  a<Jl)mT\\.  (37f) 


Except  for  the  three  exponential  degradation  factors,  the  above  estimate  is  equal  to  one  (the 
desired  value).  The  first  two  exponential  factors  degrade  the  estimate  as  a  consequence  of  sampling  the 
NCE  function  at  discrete  points  along  the  r  and  a  axis  of  the  ambiguity  plane.  This  results  in  a  two* 
dimensional  scallop  loss  (in  the  estimate)  over  the  ra -plane.  When  C|  and  «2  are  zero,  these  factors 
reduce  to  one.  The  third  exponential  factor  degrades  the  estimate  as  a  consequence  of  the  imperfect 
method  of  compensating  for  time  scale-factor  compression  (or  expansion)  in  the  algorithm.  The  design 
considerations  which  may  be  employed  to  limit  the  degradation  of  the  NCE  estimate  will  now  be 
addressed. 

Expected  Degradation  of  the  NCE  Estimate 

The  NCE  estimate  given  by  Eq.  (37)  is  a  function  of  the  variables  €,,  e2,  a0,  v,  and  A,  and  the 
design  parameters  M  and  Bj.  The  variables  e(,  e2,  and  v  are  zero-mean  random  processes,  while  A  is 
either  constant  or  is  comprised  of  a  mean  value  with  a  random  component  (over  the  analysis  interval). 
For  purpose  of  the  analysis  to  follow,  it  will  be  reasonable  to  assume  that  the  random  processes  are  all 
ergodic  and  statistically  independent. 

A  study  of  Eqs.  (36)  and  (37)  reveals  that  the  amplitude  function  A  (mT,)  can  significantly 
influence  the  NCE  estimate  if  this  function  is  highly  unstationary.  For  example,  suppose  that  A  (mT,) 
for  a  particular  m  is  much  greater  than  M  times  the  value  for  the  remainder  of  the  set  of  m.  In  this 
event,  the  value  of  the  NCE  estimate  given  by  Eq.  (36)  will  approximately  equal  one  regardless  of  the 
phase  variation  over  m.  This  is  a  nontrivial  problem  and  has  occurred  in  practice  in  connection  with 
transient  signals  of  short  duration  (comparable  to  the  sample  period  T,).  It  can  readily  be  perceived 
that  the  effect  of  relatively  high-level  transient  bursts  in  the  amplitude  level  is  to  shorten  the  effective 
(or  useful)  integration  time  of  the  NCE  estimate.  This  in  turn  reduces  (rather  than  increases)  the 
degradation  effects  under  consideration.  Consequently,  to  obtain  a  useful  measure  of  the  degradation 
effects  of  the  phase  parameters  under  consideration,  the  amplitude  parameter  Rm  will  be  considered  as 
constant  over  the  integration  time.  This  procedure  will  maximize  the  phase  misalignment  degradation 
effects  of  the  NCE  estimate,  which  is  of  primary  concern  in  this  paper.  Therefore  letting  Rm  —  1,  the 
expected  degradation  will  be  determined  for  each  of  the  three  factors  in  Eq.  (37)  separately. 

Degradation  Due  to  Sampling  Error  * 2 

Consider  first  that  <i  and  a0  are  both  zero.  The  value  of  the  NCE  estimate  due  to  the  error  c2 
becomes  (23] 

Xk(T0.-*2Bj/2M/c)  -  \j-z  I  *  “  I 

•  M  m—Mli 


sin  («jtr/2) 

M  sin  (cjir/2Af) 


sin  (ntj/2) 

w«i/2 


<  —  sin  (»/4)  -  0.90. 
w 


(38a) 
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Thus,  the  expected  degradation  due  to  imperfect  a  alignment  will  be  less  than  about  0.91  dB.  Further, 
since  the  probability  density  function  for  e2  will  be  constant  over  the  range  -1/2  <  «2  <  1/2,  the 
expected  degradation  (averaging  over  e2)  will  be 


i*i/2  sin  (irg2/ 2) 

•*®  ire2/2 


de2 


2U2J4L  =  0.96  (-0.33  dB). 
rr/4 


(38b) 


Thus,  the  increment  sample  size  along  the  a  dimension  appears  to  be  suitably  chosen  for  practical 
applications. 


Degradation  Due  to  Sampling  Error  t  ( 


Consider  next  that  «2  and  a0  are  both  zero, 
closely  approximated  as 


■  1  a 

M  **  e 

m — M/2 


Assuming  that  1  <  <  M,  the  NCE  estimate  may  be 


e 


i2nt  |*/(/ )/ Bj 


dt\ 


cos  \2irt\vUM Bj]dt  +  / 


sin  \2ire\v(tM Bj]dt\. 


(39a) 


Since  v(t )  is  a  zero-mean  function  over  the  analysis  interval,  the  imaginary  term  of  the  NCE  estimate 
will  be  near  zero  and  the  real  term  can  be  expected  to  dominate  over  the  permissible  range  of  e  t  and 
i >(t).  Although  the  expected  value  of  the  imaginary  term  will  be  zero,  it  does  not  follow  that  the 
expected  value  of  the  NCE  estimate  will  be  determined  by  the  real  term  alone.  However,  it  does  follow 
that  a  lower  bound  on  the  expected  value  of  the  estimate  can  be  determined  by  using  only  the  real 
term.  And  as  long  as  the  real  term  does  not  become  small  compared  to  one,  it  will  closely  approximate 
the  true  expected  value  of  the  estimate.  Thus, 

1  p  T/2 

X  (t,0)  >  |—  J  cos  litres (t)/Bj\dt\.  (39b) 

And  letting  p„U>)  be  the  probability  density  function  of  v  and  assuming  that  p„(—v)  —  pv(v),  the  lower 
bound  on  the  expected  value  of  the  NCE  estimate  (when  ct  is  given)  is 

£{i  (t,0)}  >  2f  p„(*)cos  (7irt\vl Bj)dv  -  2  f  p«(0)cos0  d9  (39c) 

<c0  •'O  *00 

where  0  -  2 ir« ,»»/£>  and  p„(9)  -  (Bjl2vtx)pv(Bfil2vt\). 


The  above  relation  has  been  solved  for  three  probability  density  functions,  and  the  results  are 
tabulated  in  Table  3.  The  first  probability  density  is  the  case  where  v  is  uniformly  distributed  over  a 
bandwidth  Bp.  The  second  is  a  Gaussian  function  whose  standard  deviation  is  limited  to  a  maximum  of 
Bj/4.  (The  bandwidth  Bv  is  the  information  bandwidth  of  the  signal  (see  Ref.  30,  pp  229-236).)  The 
probability  density  for  the  third  case  is  realized  when  v  is  a  sinusoidal  fluctuation  whose  peak-to-peak 
frequency  excursion  is  uniformly  distributed  over  the  bandwidth  Bv  (33).* 


•The  three  probability  density  functions  considered  here,  as  well  as  a  number  of  other  distributions,  may  be  found  in  Ref.  33 
This  earlier  work  demonstrates  that  the  coherence  degradation  will  depend  essentially  on  the  standard  deviation  of  the  #  variable 
(and  be  relatively  independent  of  the  probability  density  function)  when  it,  is  less  than  one. 
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Table  3— Probability  Density  and  Expected  NCE  Estimates 
Relating  to  the  Sampling  Error  <  | 


PROB.  DENSITY  />» 

<tJB„ 

£{x*ft(r.0)) 

1 

1  IB, 

<1*1  <  BJ2) 

1/2V3 

Bv  <  Bj 

sin  (ntiBJBj) 
irt\BJBj 

2 

1  -**2ib2 

B,  f 

(kl  «  Bj/2) 

1/V2? 

B ,  <  yfnTl  Bj 

-nlt.BJB,)2 

e  ‘  ' 

3 

'IMMmiSEBBSii 

1/2V6 

Bv  <  Bj 

Bj 

irfi,  ‘  2 1  v  | 

(kl  <  BJ2) 

*<i8,  Jo 

-f 

~e 

From  the  table  one  can  observe  that  the  expected  degradation  will  depend  on  B„  and  the  distribu¬ 
tion  of  v  over  Bj,  as  well  as  the  sampling  error  <<.  For  |«||  —  1/2  and  B„  at  its  maximum  permitted 
value,  the  lower  bound  of  the  expected  degradation  in  each  case  will  be:  0.637  (-3.92  dB)  in  case  1, 
0.735  (-2.68  dB)  in  case  2,  and  0.814  (-1.79  dB)  in  case  3.  Since  the  error  will  be  uniformly  distri¬ 
buted  over  the  range  -1/2  <  e,  <  1/2,  the  expected  degradation  due  to  this  cause  in  each  case  (averag¬ 
ing  over  c  |)  will  be  limited  to 


SijirBJlBj)  Si(ir/2) 
nBJ2Bj  *  n/2 


0.873  (-1.18  dB), 


and 


~~  erf{~-  <  2\/Fct/(47I')  “  0  905  (_0  86dB)> 

2V?  T' T}  «  - 0,3 


where  the  error  function  erf  (z)  is  defined  as  (Ref.  24,  pp.  295-300) 

erfiz)-^Soe~xtdx- 

This  amount  of  degradation  appears  tolerable  for  practical  application.  However,  the  expected 
degradation  can  be  further  reduced  by  sampling  the  one  channel  at  twice  the  Nyquist  interval  and  inter¬ 
lacing  the  sample  sets  in  Eq.  (30)  as  described  earlier.  In  this  event.  Table  3  will  still  be  applicable  with 
the  understanding  that  t|  is  limited  to  ±  1/4. 


NRL  REPORT  8438 


Degradation  Due  to  Imperfect  Doppler  Compensation 


Finally,  consider  that  and  e2  are  both  zero, 
situation  may  be  written 


M/2-1 


1  M/Z- 

Xjr0.a0)  -  I—  £ 


-ibraQmT 


«— MI2 


Assuming  1  <  <  M,  the  NCE  estimate  in  this 


1  f  r/2  rTI 2 

—  -j  !  J _r/2Cos  [2ira0/p(r)ldr  +  ij  r/2sin  [2wa0/»'(r))d/| 

.1  /*r/2 

>  ly  j_7./2cosl2ira0r»' (/)]<* I,  (40) 


where  T  *■  (A/  -  1)7",  —  ( M~l)/Bj .  Again,  letting  ^(v)  be  the  probability  density  of  v  and  assum¬ 
ing  that  pv(~f)  “  Pf(u),  the  lower  bound  on  the  expected  value  of  the  NCE  estimate  becomes 

E[xkQ(r0ta0))  ^  —  Jo  />„(*)!  J0  cos(2ira0fi/)<//|</i/ 


-  JitVwi 


sin  (iraoTv) 


ira^Tv 


dv 


>  2 /;„» 


sin  (ira^Tv) 
ITIXqTv 


dv 


/>„(*<)  cos  (iraoTvxOdvldlx. 


(41) 


One  may  recognize  that  the  inner  integral  in  Eq.  (41)  is  identical  in  form  to  Eq.  (39c).  Conse¬ 
quently,  the  results  given  in  Table  3  will  be  applicable  to  this  case  provided  a0B„Tx/2  is  substituted  for 
e \BJBj.  The  lower  bound  on  the  expected  values  of  the  NCE  estimate  (for  the  three  given  probability 
densities)  may  therefore  be  computed  as 


SHira0B,T/2) 
na0B,T/2  * 

(42a) 

erfiJZaoB'T/2) 

<*o BJ 

(42b) 

erf(wa0Bv  T/4VJ) 
y/ira0B,T/2yn  ' 

(42c) 

Since  the  above  expectations  are  derived  from  ensemble  averages  for  the  random  variable  v(r),  it 
will  be  informative  to  determine  the  error  in  the  NCE  estimate  for  specific  examples  of  v(t)  for  com¬ 
parison  purposes.  Two  examples  are  chosen  which  produce  rather  severe  degradation  on  the  estimate. 
In  the  first  example  let 

vit)  -  ±  B¥t/T  (-T/2  <  /  <  r/2)  (43a) 
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and  in  the  second  example  let 


At)- 


±  BJ2  (for  0  <  /  <  772) 

*  BJ2  (for-T/2  <  /  <  0)  • 


In  the  first  example,  the  frequency  varies  linearly  from  7  BJ 2  to  ±  BJ 2  over  the  integration  period  T. 
This  is  representative  of  what  can  occur  in  practice.  However,  it  is  a  severe  example  in  that  v(t)  is  per¬ 
fectly  correlated  with  the  integration  variable  t  (which  will  maximize  the  degradation  due  to  the  linear 
frequency  slide).  In  the  second  example,  the  frequency  v(t)  remains  fixed  at  ±  BJ 2  for  one-half  the 
integration  time  and  then  flips  to  the  negative  value  for  the  remainder  of  the  integration  period.  This  is 
an  extreme  case  which  will  be  approached  infrequently  in  practice.  However,  it  will  provide  a  suitable 
upper  bound  on  the  expected  NCE  degradation. 

Employing  the  two  specific  examples  in  Eq.  (40)  and  carrying  out  the  integration  gives 

-  »  c(VOTT) 


X.  (t0,oo)  - 


for  the  first  example,  where  C(— )  is  the  Fresnel  cosine  integral  of  the  indicated  argument  (Ref.  24,  pp. 
300-304),  and 


X*  (r0,ao) 
*0 


sin  (tra0ByT/2) 
naoBvT/2 


for  the  second  example.  Graphs  of  these  functions  along  with  those  given  in  Eq.  (42)  are  shown  in 
Fig.  6. 

The  curves  (Fig.  6)  illustrate  that  the  value  of  a0B„T should  be  limited  to  about  one  if  excessive 
degradation  due  to  imperfect  Doppler  compensation  is  to  be  avoided.  Since  all  of  the  curves  are  above 
•1  dB  for  a0B„T  equal  to  (or  less)  than  0.5,  this  value  would  represent  a  conservative  choice  for 
a0B,T  .*  Therefore,  a  reasonable  upper  bound  on  the  integration  parameter  M is  (from  Eq.  (31)) 

^<-^-^-  +  1.  (45a) 

l«*0l  B* 

And  since  the  Doppler  ratio  a0  is  Avjc,  where  &v0  is  the  difference  in  source  speed  along  the  propaga¬ 
tion  paths  to  the  two  signal  sensors  and  c  is  the  signal  propagation  speed  in  the  transmission  medium 
[30], 

+  +  <45b, 

In  the  case  of  underwater  acoustic  applications,  c  is  approximately  2880  knots.  Assuming  a  source- 
speed  differential  of  10  knots,  a  suitable  bound  for  M  is  less  than  or  equal  to  288  BJB,. 

SUMMARY  AND  CONCLUSIONS 

Although  the  Fourier  Transform  of  a  temporal  function  is  normally  used  to  decompose  the  func¬ 
tion  into  its  complex  spectral  components,  the  sectionalized  Fourier  Transform  (SFT)  may  be  employed 

•W.H.  Marsh  (191  has  recommended  that  atB,T  be  limited  to  leas  than  0.2$  in  connection  with  the  MSC  estimate.  However  thia 
appears  to  be  ultra-conservative  in  lisht  of  the  results  displayed  in  Fig.  6. 
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Fig.  6  —  Coherence  degradation  due  to  imperfect  Doppler  compensation 
as  a  function  of  a0ByT.  (Curves  labeled  I.  2.  and  2  result  from  assuming 
the  corresponding  probability  densities  given  in  Table  3) 


in  a  manner  which  preserves  the  temporal  properties  of  the  original  signal.  It  has  been  demonstrated 
that  when  a  contiguous  sequence  of  J  SFT  spectral  bins  is  appropriately  summed  (J  typically  being  a 
small  number),  the  resulting  transform  of  a  bandlimited  signal  is  equivalent  to  basebanding,  filtering, 
and  sampling  the  signal  in  the  time  domain.  The  value  of  /  may  be  chosen  to  control  the  passband  and 
leakage  characteristics  of  the  filter  (see  Figs.  1-4). 

With  the  advent  of  the  FFT  algorithm  and  modem  array  processors,  use  of  the  SFT  to  baseband, 
filter,  and  sample  signals  considerably  simplifies  the  programming  of  multidimensional  correlation  pro¬ 
cessors  in  practical  applications.  Further,  using  the  ambiguity  kernel  as  an  approximation  of  signal  time 
compression  (or  expansion),  the  FFT  algorithm  is  applicable  to  correlation  mapping  along  the  Doppler- 
ratio  axis  of  the  ambiguity  plane.  The  resulting  error  has  been  shown  to  be  negligible  when  the  product 
of  the  signal  bandwidth  and  the  correlator  integration  time  is  less  than  the  inverse  of  the  maximum 
Doppler  ratio  being  employed.  Using  the  techniques  described  in  this  paper  (and  modification  thereof), 
two-dimensional  correlation  mapping  of  low-frequency  acoustical  signals  over  long  integration  intervals 
has  been  implemented  for  the  NCE  estimate  (or  the  MSC  estimate)  well  in  excess  of  real  time. 

It  may  be  concluded  that  the  sectionalized  Fourier  Transform  has  many  applications  as  an  alter¬ 
nate  (and  convenient)  method  of  time-domain  processing  using  modern  array  digital  computers.  Its 
application  is  limited  only  by  the  sample  rate  which  can  be  processed  in  the  digital  computer  employed. 
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Appendix 

EXPECTED  VALUE  OF  THE  DISTORTION  FUNCTION 


From  Eqs.  (3)  and  (5),  the  distortion  function  may  be  written  in  the  form 


where 


and 


J>  -  jr  C’, 

(Ala) 

Am  -  A,  (mf|  +  t)/A)(mT\), 

(Alb) 

l  r1 

Vm  “  7  Jo  »kl"iTt  +  x)dx, 

(Ale) 

F(t)  —  sin  (wJt/Tj)/ sin  (irf/r2), 

(Aid) 

and  with  the  additional  provisions  that  A,  (mT\)  >  0  and  ±  1/2  is  added  to  8r  when  /  is  even  (as 
described  in  the  text  following  Eq.  (5». 

The  function  F(t)  is  deterministic,  while  the  functions  Am  and  v„  are  comprised  of  sample  func¬ 
tions  ( A,  and  vk)  drawn  from  ensembles  (or  sets)  with  common  statistical  characteristics.  To 
emphasize  this  fact,  Eq.  (Ala)  is  rewritten  as 

D.,  -  (A2> 

where  the  subscripts  j  and  k  have  been  attached  to  the  sample  function  A„  and  vm  to  imply  that  a  sam¬ 
ple  member  from  each  ensemble  M,}  and  |v  J  is  chosen  in  the  computation  of  the  distortion  function. 

The  expected  value  of  the  above  distortion  function  is  obtained  by  computing  the  integral  aver¬ 
aged  over  the  double  ensemble  of  sample  functions.  In  computing  the  ensemble  averages  it  will  be 
assumed  that  the  two  sample  sets  are  independent,  and  the  sample  functions  At  and  vk  from  each 
ensemble  are  both  ergodic  and  stationary  over  the  integration  interval  T2.  (The  condition  of  stationary 
does  not  apply  to  the  functions  Ami  and  pmk  which  are  constructed  from  At  and  vk.)  Restrictions  on 
the  properties  of  A,  and  v*  to  qualify  as  a  member  of  each  ensemble  will  be  determined  as  we  proceed 
with  the  analysis.  (Basically,  the  members  of  each  ensemble  must  be  such  that  their  interrelated  power 
spectral  density  is  confined  to  fall  essentially  within  the  spectral  window  of  J/T2  Hz.)  Under  the  stipu¬ 
lated  conditions,  the  expected  value  of  the  distortion  function  becomes 

E{DJ  -  <DmJk>  (averaged  over  j,k) 

-  ~r  f  <  AmJ>  Fit)el2w>(,/Ti<  e2ui’*'>  dt  (A3a) 

#2  §r* 

Letting  P\(Am,t )  and  p2(vm,t)  denote  the  time-variable  probability  densities  of  Am  and  vm  respectively, 
the  respective  ensemble  averages  can  be  written 

<4w>  -  /0  AmPMm.t)dAm  (A3b) 
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and 

<ell'"'mk'>  ~  $_"Pi Gm-i)*'2w*m'di>m.  (A3c) 

At  this  point  it  is  important  to  realize  that  the  restrictions  on  the  member  functions  of  the  two 
independent  ensembles  are  interrelated.  The  functions  [Am,]  are  amplitude  modulations  of  the  carrier 
functions  exp(/2ir(8(./T2  +  »'**Jfl-  A  member  of  this  carrier  set  is  exp(<2ir[6f/7’j  + !>*#]/),  where  vu  is 
the  extreme  limit  of  a  slowly  varying  fluctuation  frequency.  Therefore,  since  any  dynamic  variation  in 
the  sample  function  Amj  will  produce  equal  sidebands  of  power  on  either  side  of  the  carrier  frequency 
bjT2  +  I'm,*  the  carrier  frequency  must  be  restricted  to  lie  sufficiently  well  within  the  spectral  window 
to  accomodate  the  amplitude  modulation  sidebands.  On  this  basis  it  is  evident  that  when  J  is  one  and 
8,.  —  1/2,  Am  must  be  essentially  constant  over  Tj  if  the  distortion  function  is  not  to  become  exces¬ 
sive.  On  the  other  hand,  when  J  is  large,  the  parameter  8(.  will  not  play  a  significant  role  in  restricting 
the  members  of  the  ensembles.  Consequently  for  1  <  <  J,  the  spectral  window  need  be  only 
sufficiently  broad  to  accomodate  both  the  maximum  deviation  in  the  frequency  fluctuation  v(t)  and  the 
spectral  sidebands  introduced  as  a  consequence  of  the  amplitude  modulation  (if  significant  distortion  of 
the  resulting  transformed  signal  is  to  be  avoided). 

Since  the  two  ensembles  of  sample  functions  are  assumed  independent,  it  will  be  convenient  and 
appropriate  to  treat  the  amplitude  and  frequency  fluctuation  problems  separately. 

Expected  Distortion  Due  to  Frequency  Fluctuations 

Consider  first  that  the  functions  [Aj(mT\  +  f)}  are  essentially  constant  over  time  increments  of  T2 
seconds  such  that  the  ensemble  [Ami]  is  unity.  Considering  then  the  broad  class  of  ensemble  functions 
vmt,  it  is  reasonable  to  assume  that  for  every  member  function  there  exists  a  member  function  which 
has  its  negative  time  characteristics.  This  is  equivalent  to  assuming  that  p2(vm.  - 1)  -*  p2(v„,i).  With 
these  considerations,  Eq.  (Ala)  reduces  to 

£|0J  -  -jr  J*0  ^  gjn C0S,2w(6c  +  T^m)t/T2\dvndi  (A4) 

It  is  well  to  note  that  the  maximum  time  we  need  be  concerned  with  is  T-Jl.  However,  the  mag¬ 
nitude  of  the  Dirichlet  kernel  decays  rapidly  for  /  greater  than  T-J1J  —  7'l/2.  Consequently,  the 
significance  of  p2{Vm,t)  becomes  increasingly  less  beyond  t  ~  T,/2. 

Since  vmk  is  the  running-time  average  of  i'/k(m7'|  +  r),  its  probability  density  function  wil!  depend 
on  the  dynamic  characteristics  (or  power  spectrum)  of  v(i)  as  well  as  on  the  smoothing  time  t.  It 
should  also  be  apparent  that  the  peak  magnitude  of  v„k  cannot  exceed  the  peak  magnitude  of 
vk(mT\  + 1)  (that  is,  lvmkl  <  lvk(mT\  +  r)|).  When  t  is  sufficiently  small  so  that  vk(mT\  +  /)  is 
approximately  equal  to  vk(mT\)  +  vk(mT\)t  over  all  members  of  the  ensemble,  the  probability  density 
of  vm  will  approximate  the  probability  density  of  v.  On  the  other  hand,  when  t  is  large,  the  probability 
density  of  »m  will  be  compressed  relative  to  the  probability  density  of  v.  This  is  evident  since  the  aver¬ 
age  of  a  rapidly  varying  zero-mean  function  approaches  zero  over  relatively  long  time  intervals. 

To  exemplify  the  above  consideration,  let 

pk(mT\  +  t)  ■■  vwsin  (2irpt  +  dm) 


*M.  Schwartz,  W.R.  Bennett,  and  S.  Stein,  Cammunicatiom  Svitems  and  Tcrlmiqurs,  McGraw-Hill  Book  Co..  N  Y..  I%6.  pp  173- 
176. 
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where  vM  is  the  frequency  deviation  <  fly/2  -  J/2T7)  and  p  is  the  modulating  frequency.  Then 

-  “X  *'*<*7'i  +  *><**  “  •'w  S'npt"  sin  (*p'  +*») 


-  £-"  ”pt  vk(mT\  +  /  -  t/2). 
irpt 

Thus,  vmk  is  a  delayed  (by  the  amount  t/2)  and  compressed  version  of  the  function  vk(mT\  +  t).  The 
ratio  vu/p  is  known  as  the  frequency  modulation  index  of  the  process,  and  the  resulting  signal  spectral 
energy  can  be  expressed  in  terms  of  Bessel  functions.*  To  prevent  significant  spectral  energy  of  the  sig¬ 
nal  from  exceeding  the  passband  of  the  spectral  window,  this  ratio  should  be  no  greater  than  n/2.f 
This  places  an  upper  bound  on  the  modulating  frequency  p  of  2vuht  <  Bj/v. 


Using  the  upper  bound  for  p,  the  function  vmk  becomes 


sin  (2 vstt) 
2 vMt 


vk(mT\  +  t  -  t/2). 


and  for  |/|  less  than  or  equal  to  7V2,  the  compression  factor  is 


sin  (2i lyt) 
2  vm!) 


sin  vu,T i 

*>uT\ 


<  2  sin  (1/2)  =  0.96. 


Thus,  in  the  extreme  case  (for  p  maximum  allowable),  the  probability  density  for  vm  will  be  approxi¬ 
mately  the  same  as  the  probability  density  for  v,  within  the  primary  lobe  of  the  Dirichlet  kernel.  The 
maximum  rate  of  change  of  vk(mT{  +  t)  is  2irpvM  <  Av\,  <  B}  -  fi/T}.  This  rate  limits  the  change 
in  vk(mTy  +  t)  over  the  period  of  Tt  seconds  to  less  than  the  spectral  window  width  of  fly  Hertz  (or 
equivalently,  vT\  <  1).  Little  error  will  therefore  result  in  assuming  that  the  probability  density  of  vm 
is  the  same  as  the  probability  density  of  v  over  T2. 

As  a  consequence  of  the  above  analysis,  Eq.  (A4)  will  closely  approximate  (using  Eq.  (6a)) 

£{D-1  =  1 cosln(6‘  +  TphlMr 


“  Wj(&c  +  TiV)dv, 


(A5) 


where  p„(k)  is  the  probability  density  function  for  v  and  Wj(-)  is  the  spectral  window  function  for  the 
indicated  argument. 

The  expected  value  of  the  distortion  function  is  therefore  the  weighted  value  of  the  magnitude  of 
the  spectral  window  function  (see  Figs.  1  and  2  of  main  text).  The  weighting  function. is  the  probability 
density  function  for  the  fluctuating  frequency  v.  When  J  >  >  1,  Wj(x)  is  approximately  equal  to  one 
over  the  spectral  window  passband.  Therefore  if  p„(p)  is  zero  for  Ik  I  >  fly/2,  the  expected  value  of 
the  distortion  function  is  essentially  unity. 


’op.  cit.,  pp  225-221 
fop.  cH„  pp.  120-121 
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Expected  Distortion  Due  to  Amplitude  Fluctustions 

To  study  the  expected  distortion  due  to  amplitude  modulation  we  shall  consider  that  v(mT\  +  i) 
is  constant  over  -T-jl  <  t  <  7V2,  and  let  xm  —  8r  -4-  T2v(mTx).  The  ensemble  average  will  be 
taken  only  over  j  holding  k  fixed,  so  that  Eq.  ( A3a)  becomes 

£*(£>„}  “  < Dmjk>  (averaged  over  j) 


J_  CT ^ 

r2 


<*mi> 


sin  (nJt/T2)  /2„, t/T,  . 

- e  at. 

sin  (ir//r2) 


(A6) 


To  proceed,  it  will  prove  convenient  to  perform  the  ensemble  averaging  in  steps.  First,  it  may  be 
assumed  that  for  every  member  function  Ami  there  exists  a  complement  member  function  with  nega¬ 
tive  time  characteristics.  Thus,  the  average  of  the  member  and  its  complement  is  an  even  function  of 
time.  As  a  consequence,  the  expected  value  of  the  distortion  will  be  real  and  we  need  carry  the  averag¬ 
ing  process  only  over  even  functions  of  time.  The  symbology  Amt  will  henceforth  be  used  to  represent 
an  even  function  of  time,  with  the  understanding  that  the  first  step  of  ensemble  averaging  has  been 
effected. 


Next,  let  At(mT\  +  i)  be  written  as  /f0,U  +  a,(/)],  where  A0/  is  the  mean  value  over  the  time 
interval  T2  and  a,  is  a  zero-mean  function  greater  than  -1  over  this  interval.  Then  Ami  takes  the  form 
(1  +  <2/0)1/11  +  0/(0)).  Over  the  time  interval  7j,  o, (r)  may  be  expanded  into  the  Fourier  series 

p 

a j{t)  -  £  b,/,  cos  (2 irpi/T2) 

where 


<*,«>)-  £  bip 


p-i 


>  -1. 


Since  no  significant  spectral  energy  will  be  permitted  to  fall  outside  of  the  spectral  window,  the 
ensemble  of  functions  a,  to  be  considered  wilt  be  limited  to  those  whose  upper  limit  P  is  restricted  by 
the  relation 

P^J/2-  UM| 

where  Uw|  is  the  maximum  excursion  of  the  frequency  fluctuation  xm.  This  informs  us  that  when  7  is 
either  less  than  2  or  |xw|  is  7/2,  no  significant  amplitude  modulation  can  be  permitted  without  serious 
distortion  of  the  resulting  transformed  signal.  If  one-half  of  the  spectral  window  is  reserved  for  fre¬ 
quency  modulation,  /*  will  be  limited  to  values  less  than  7/4. 


With  the  above  considerations  then,  Eq.  (A6)  reduces  to 


£*U>J  -  <  J0' 


l  +  £  bip  cos  (npt) 


1  +  a,(0) 


sin  (ir7r/2) 
sin  (ir//2) 


cos  (irxmt)dt  > 


“  <  l  £  */|,i  VAxm  +  p)V  £  b, |,|  > .  (A7) 

p—p  p—p 

where  bl0  ”  2  and  W ](•■■)  is  the  spectral  window  function  (defined  in  Eq.  (6))  for  the  indicated  argu¬ 
ment. 


The  above  relation  shows  that  the  expected  signal  distortion  due  to  amplitude  fluctuations  will  be 
dependent  on  the  flatness  of  the  spectral  window  over  the  passband.  Thus,  when  7  is  sufficiently  large 
so  that  the  spectral  window  can  be  assumed  to  be  unity  over  the  passband,  Eq.  (A7)  reduces  to  one. 
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